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A NEW FORMULA FOR SOME LINEAR STOCHASTIC 
EQUATIONS WITH APPLICATIONS 

By Offer Kella 1 and Marc Yor 
Hebrew University of Jerusalem and Universite Pierre et Marie Curie 

We give a representation of the solution for a stochastic linear 
equation of the form X t = Y t + f, Q ( , X s _ dZ s where Z is a cadlag 
semimartingale and Y is a cadlag adapted process with bounded vari- 
ation on finite intervals. As an application we study the case where 
Y and — Z are nondecreasing, jointly have stationary increments and 
the jumps of — Z are bounded by 1. Special cases of this process are 
shot-noise processes, growth collapse (additive increase, multiplica- 
tive decrease) processes and clearing processes. When Y and Z are, 
in addition, independent Levy processes, the resulting X is called a 
generalized Ornstein-Uhlenbeck process. 



1. Introduction. In this paper we show that when Z is a cadlag adapted 
semimartingale and Y is cadlag adapted and with bounded variation on 
compact intervals, then the unique cadlag adapted solution of X t = Y t + 

t , X s _ dZ s is given via the representation X t = Jj t i U Uj t dY u where U u j 
is defined by formula (2) below. This form seems to be new and we note 
that the integral with respect to Y is defined path-wise while the integral 
in the integral equation can be a stochastic integral. Of course when Y is 
a semimartingale, one cannot expect such a representation of the solution 
since {U u j\0 < u < t} is not adapted as a process indexed by u. 

We discuss an application to the case where Y and —Z are nondecreasing 
processes jointly having stationary increments and subsequently specialize to 
cases where one or both also have independent increments (Levy processes). 
This model is a generalization of both the shot-noise process as well as 
a growth-collapse process (e.g., see, [7, 11, 16] and references therein) or 
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more generally an additive increase and multiplicative decrease process. The 
later have been used as models for the TCP window size in communication 
networks. 

We note that Jacod ([8], Theorem 6.8, page 194) and Yoeurp and Yor 
[21] give a complete solution for the case where the integrator is a semi- 
martingale and the driving process is cadlag, Jaschke [9] gives a derivation 
for the case where the integrator does not have jumps of size —1, and Protter 
([20], Theorems 52 and 53, pages 322-323) treats the case with a continuous 
integrator. 

The literature related to generalized Ornstein-Uhlenbeck processes and 
their applications which are directly related to some of the special cases of 
the applications that we consider is huge and growing exponentially fast. We 
refer the reader to [1-6, 14, 15, 17-19, 22] and further references therein. 

2. Main result. With respect to some standard (right continuous aug- 
mented) filtration, let Y = {Yt\t > 0} and Z = {Zt\t > 0} be two adapted 
cadlag processes. Denote Zq^ = 0, and for t > 0, Z t ~ = lim s -^ Z s . Set AZ t = 
Zt — Zt— when Z is of bounded variation on compact intervals (BV); set 
Z$ = Zt — Xlo<s<t and similarly for any other cadlag process consid- 
ered in this paper. 

Theorem 1. Assume Y and Z are cadlag and adapted, Y is BV and Z 
is a semimartingale. Then the unique cadlag adapted solution to the equation 
X t = Y t + J (0jt] X s _ dZ s is 

(1) X t = I U u>t dY u , 

J[o,t] 

where 

0<u<t, 
0<u = t 

and [Z,Z] is the quadratic variation process associated with Z. When Z is 
BV then (2) reduces to 

( e z c t-zz ■ [ (i + AZ s ), 0<u<t, 

( 3 ) U U,t = < u <s<t 

[l, 0<u = t, 

where Z c is the continuous part of Z as defined earlier (rather than the 
continuous martingale part of Z as is customary in stochastic calculus). 



(2) 



U u , t 



[ (l + AZ s )e~ 



AZ S 



u<s<t 



u, 
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Proof. Note that with T = and for n > 1, T n = inf{t > T n _ x \AZ t = 
— 1}, then for T n < u < t < T n+ \ 

(4) ^^ = U Ujt (l + AZ u ), ^£ = U Ujt . 

Also, since Y is a BV process, the covariation process [1", Z] is given via 
[Y,Z] t = ^2 0<s<t AY s AZ s . If one follows the solution in equation (6.9) in 
Theorem (6.8) on page 194 of [8], then for T n <t < T n+ \ we have that 



x t = u Tntt Ur Tn + I Uj* dY u - [ U^ u d[Y, Z] u ) 

\ J(T n ,t] J(T n ,t] J 

(5) =U Tn . t AY Tn + f U u , t (l + AZ u )dY u - U Uit AY u AZ u 

J(T„,t] 



'{Tn,t] J(T n ,t] 

'( T "'*l T„<u<t 
U u ,t dY u , 

[T n ,t] 

where the second equality is justified since the first integral on the right-hand 
side of the first equality is a path- wise Stieltjes integral, and the second is a 
sum which is also defined path-wise. If Y was a general semimartingale, then 
interchanging Ux ni t with the integral sign like this would not be justified as 
the resulting integrand would no longer be adapted. Clearly if n > 1, then 
U u ,t = for u < T n , and thus 



(6) X t = / U U)t dY u . 

J [OA 

Since this holds for all n, the proof for the more general case is complete. 
For the case where Z is BV, it is evident that [Z, Z] c = 0, and it is easy to 
check that ^2 u<s<t AZ s is convergent (actually, absolutely convergent), and 
hence the result follows. □ 

Of course one may also define the counting process, 

(7) #t= E k^-i}- 

0<s<t 

which is a.s. finite for all t > and right-continuous (possibly a.s. identically 
zero or terminating), and write 



(8) X t = / U u , t dY u . 

J[ T N t A 

It is worth while to note that for the case where Z is also a BV process, 
there is a more direct proof involving (path-wise) Stieltjes integration which 
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can be taught in a classroom as follows. Write Z = A — B, where A and B are 
right-continuous and nondecreasing and have no jump points in common. 
Write A^ = At — A1 = X^o<s<t max (AZ s , 0) and similarly for B. Observe that 
by right continuity A At, ABt, Af — Aq and Bf — Bq all converge to zero as 
t \. 0. In particular, for every t for which —1 < AB S (< 0) for < s < t, we 
have that 



(9) 1 + A?-A < ' [ (l + AA s )<e 



0<s<t 

and 

(10) 1 + Bf - B < + AB S ) < e B ?- Bo 

0<s<t 

which implies that 

(11) n (i + az s ) = ( n a + AA a )\ ( n (i + ^b s 

0<s<t V)<s<t ' V)<s<i 

as 1 10. 

Now note that with Ct = e z t and Dt = Y\_Q <s<t (\ + AZ a ) , ordinary (Stielt- 
jes) integration by parts yields 

U t = C t D t = C 0+ D 0+ + [ D s _dC s + [ C s .dD s + V AC S AD S , 

J(P,t] J{o,t] 0<s < t 

(12) 

and it is easy to check that the continuity of C and the fact that dCt = Ct dZ^ 
imply that 



(13) U t = l+ / U a -dZ s . 

J(o,t] 

With this formula established, it is clear that if we denote U u t as in (3), 
then in an identical way to which (13) was obtained we have (path-wise) 
that 



(14) U U)t = l+ / U u ^dZ s 

J(u,t] 

for all < u < t. 

Now, if X t = Jj 0jt , U S)t dY s , then X t _ = Jj Q ^ U s>t - dY s and thus X s _ dZ s 
is given by 



U U)S - dY u dZ s = / U Uj8 - dZ s dY u 

(0,t]J[0,s) J[0,t)J(u,t] 



(15) 

i[o,t) 



[ (U U:t -l)dY u , 
J\o,t) 
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but since Ut t = 1 we can include t in the domain of integration without 
changing the value which gives 

(16) f X s _dZ s = I (U Uit -l)dY u = X t -Y t 
J(o,t] J[o,t] 

as required. 

3. Applications. Assume that Y and Z are right-continuous and nonde- 
creasing jointly having stationary increments in the strong sense that the 
law of 8 S ( Y, Z) is independent of s where 

(17) 6 a (Y(t), Z(t)) = (Y(t + s)- Y(s), Z(t + a) - Z(s)). 

It is standard to (uniquely) extend (Y, Z) to be a double sided process having 
stationary increments, that is, that t £ R rather than t > 0, thus we assume 
it at the outset. Finally we assume that Z has jumps bounded by 1. Without 
loss of generality let us assume that Yq = Zq = 0, otherwise we perform what 
follows for Y — Yq and Z — Zq which also have stationary increments. We 
consider the unique process X defined via X t = Xq + Y t — f, Q ^ A s _ dZ s for 
t > where Xq is almost surely finite; the unique solution of which is 

(18) X t = X e~ z ? II (1-AZ S )+ f e -W- z $ ]J (l-AZ s )dY u , 

0<s<t ^(°' t ] u<s<t 

where an empty product (when u = t or when t = on the right) is defined 
to be 1. 

Special cases of such processes are the shot-noise processes in which Zt = 
rt and Y are compound Poisson, growth collapse or additive increase multi- 
plicative decrease (AIMD) processes in which Y± = rt and usually Z = qN\ 
where N\ is a Poisson process with rate A, and < q < 1, as well as clearing 
processes where Z is a Poisson process or, more generally, a renewal counting 
process (see, e.g., [10, 12]). 

Consider the nondecreasing processes 

(19) Jt = Zf- Ml " AZ s )l {AZs<1} , 

0<s<t 

and Nt = Y2o<s<t ^{AZ s =i}- Then it is clear that Y, J, jointly have station- 
ary increments ^in the strong sense), and from (18) we have 

(20) X t = XQe- J n {Nt=0} + [ e-( J t~ J °h {Nt _ Ns=0} dY s . 

If J/qq q| e Js dY s is a.s. finite (recalling that for s < 0, J s < Jq = 0), then 
setting XI = fr^ t i e~( Jt ~ Ja ^{N t -N s =o} dY s it is clear that X* is a station- 
ary process. Moreover, if, in addition, either lim^oo Nt > 1 a.s. (equiva- 
lently, T\ = inf{i|AZf = 1} is a.s. finite) or Jt — > oo a.s. as i ^ oo, then 
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\X£ — Xt\ a.s. as t — > oo, and thus for any a.s. finite initial Xq, a limit- 
ing distribution exists which is distributed like Xq. 

In fact, when Xq is independent of (Y,Z), then shifting by —t, noting 
that 9-t J s = J s -t — J-t (so that 6-tJt = 0) and similarly for N and Y, it is 
clear that Xt has the same distribution as 



(21) 



J(o,t] 



X e J - t l {N _ t=0} + f e J n {Ns=0} dY s . 
J(-t,o] 



l(~t,0] 

In particular, this implies that when Xq = 0, then Xt is stochastically in- 
creasing in t > 0. 

Let us summarize our findings as follows. 

Theorem 2. // f,^ i e Ja dY s < oo a.s., and either T\ < oo a.s. or J t — > 
oo a.s. as t — > oo, then X has the unique stationary version 

(22) X* t = [ e-( J *- J °h {Nt _ Ns=Q} dY s , 

J(-oo,t] 

and for every initial a.s. finite Xq, Xt converges in distribution to Xq. More- 
over, when Xq = a.s., then Xt is stochastically increasing in t > 0. 

We note that when (Y, Z) also have independent increments so that they 
form a Levy process, then the negative of the time reversed process is a 
left-continuous version of the forward process, and thus in this case [when 
Xq is independent of (Y,Z)], X t is also distributed like 

(23) X e- J n {Nt=oy + [ e~ J n {Ns=0} dY s 

J(Q,t] 

which is also the consequence of the usual time reversal argument for Levy 
processes. In what follows we will consider special cases of this structure. 

We observe that in the general case N is a simple (i.e., a.s. ANt € {0, 1} 
for all t) counting process associated with a time stationary point process. 
Special cases of such processes are Poisson processes and delayed renewal 
processes where the delay has the stationary excess lifetime distribution as- 
sociated with the subsequent i.i.d. inter-renewal times. We will consider this 
special case a bit later. 

3.1. EXt for independent Xq, Y , Z . Since Y has stationary increments, 
it follows that EYt = EY\t. From (21) we have that when EY\ and EXq are 
finite, then for t > 0, 

r o 

(24) EX t = EXoEe J -n {N _ t=0} + EY 1 J Ee J n {Ng=0} ds, 
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and since for s < 0, we have that J s = — (Jo — J s ) is distributed like — J_ s = 
— (J_ s — Jo), and similarly for N, we have that 

(25) EX t = EX Ee- Jt l {Nt=0} + EY 1 [ Ee~ J n {Ns=0} da. 

Jo 

3.2. EXt for independent Xq, Y , Z with Levy Z . Here Z is a subordina- 
tor with Laplace-Stieltjes exponent —n z (a) = \ogEe~ aZl where, for a > 0, 

(26) Vz (a) = c z a + [ (1 - e~ ax )u z (dx) 

7(0,1] 

with c z > and J^ j, xzv 2 (dx) < oo. Since the jumps of Z are bounded above 
by 1, then v z ({l, oo)) = 0. 

In this case Z% = c z t, N is a Poisson process with rate X = v z {l} which is 
independent of the subordinator, 

(27) Jt = c z t- £ log(l-AZ s )l {AZs<1} ; 

0<s<t 

the Levy measure of which, call it za-, is defined via i>j((a,b]) = v z ((l — 
e _a , 1 — e~ fe ]) for < a < 6 < oo and with exponent 

773(a) =c z a+ (1 - e" ax )zA,-(dx) 

7(0,oo) 

(28) 



c z a+ / (l-(l-x) a )v z (dx), 

'(0,1) 



so that for a > 0, 
(29) 773(a) + A = c z a + / (1 - (1 - x) a )u z {dx). 

7(0,1] 



We note that 
(30) 



/ min(x, (dx) = / min(— log(l — x), l)z/ z (dx), 
7(o,cx>) V(o,i) 

and since — log(l — x) < < xe for < x < 1 — e _1 , the right-hand side 
is dominated above by e J, Q jx xi/ z (dx) < 00, so that Uj is indeed the proper 
Levy measure of a subordinator. Now, for this case, Ee~ Ja = e~^i^ l > s where 

773(1) = c z + [ (1 - (1 - x)> z (dx) 
7(0,1) 

(31) 

= c z + I xu z (dx) = 77^.(0) - A 
7(o,i) 
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recalling A = u z {l}. Therefore, Ee- J n {Na=0} = e -(»?i(o)-A)s e -A S = e - v ' z (o)s 
so that in this case, since rj' z (0) = c z + J^ Q ^ xu z (dx) = EZi, (25) becomes 

(32) EX t = EX e- EZ ^ + |^(1 - e~ EZ «). 

hi Zj\ 

Recall that here Y need not have independent increments. 



3.3. Independent Xq, Y , Z with Levy Y . Since for every = to < t\ < 
■ ■ ■ <t n = t the independence between Y and Z and hence the independence 
of Y and J, yield 



E 



(33) 



exp -aX>-^-xl {JVti _ i=0} (^ -Y ti _,] 



i=i 



Z 



~[exp(-r] y (ae Jt *-i l{Ar ti _ 1= o})(*i - U-i))- 



i=l 



It thus follows, as in equation (5.9) of [13] for the more general multivariate 
case and in Proposition 1 of [19] for the case where Y and Z are compound 
Poisson, that 



E 



(34) 



exp 



e~ Js l 



{7V s =o} dF s 



This implies, as in Theorem 5.1 of [13], that the conditional distribution of 
f(o t] e ~ Ja ^-{N s =o} d^s given Z is infinitely divisible, as on the right-hand side, 
—rjy/n is also a Laplace-Stieltjes exponent of a subordinator. 

Equation (34), with £o( a ) = Ee~ aX °, a A 6 = min(a, 6), and recalling 



(35) 
yields 



Ti = inf {t| AZ t = 1} = inf {t|iV t > 0} 



Ee 



-aX t 



(36) 



E£ (ae Jt l {Nt=0} ) exp^- ^ 
^o(ae- Jt )exp^-^ 



r] y {ae 



)1 



{iV s =0} 



ds 



r^lae Js )ds Jl{Ti>t} 



+ E'exp 



ae 



ids 1 



L m<t}- 
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Clearly, when either T\ < oo a.s. or Jt — > oo a.s. as t — > oo, then 
(37) lim Ee~ aXt = Eexp(- [ * ?? y (ae" Js )ds 



o 



We now observe that if N and J are independent, as for instance in the case 
where Z is a subordinator, and iV is the counting process associated with a 
time stationary version of a renewal process the latter having inter-renewal 
time distribution F having a finite mean /i, then it is well known that ./V 
is a delayed renewal process in which the times between the (i — l)th and 
ith jumps are distributed F for i > 2 and the time until the first jump (i.e., 
the delay) has a distribution with density f e (t) = (1 — F(t))/fjb. Therefore, 
in this case, 

Eexpf- J rjy(ae~ Js )ds^ = J Eexpf - J r] y (ae~ Js ) ds^j f e (t) dt. 

(38) 

Differentiating the right-hand side of the first equality in (36) once and 
setting a = gives (25) as expected, while for the case where Xq = a.s., 
differentiating twice and setting a = yields 

(39) EX? = (r,' y (0)) 2 E(j\- J °l {Ns=0} ds} -^E j\~ 2J n {Ns=0} ds. 

3.4. EX? for independent Y, Z with Levy Y,Z and Xq = 0. We note 
that for every j3 > 0, E f* e~P J " 1{jv s =o} ds = 1 ~ e ^)+ A + - , where A = v z {l}. 
Also, note that since N u < N s for u < s, 



(40) 



e~ Js lfM„=n\ ds | =2 / / e" Js ~ J »l, v=n\ duds 



{jv s =o}Qs =z / / e M^^o} 
/ Jo Jo 



e 



" (Js " Ju) e- 2J n {JVs=0} duds 



'0 Jo 

and therefore (using Fubini and the stationary independent increments prop- 
erty of J), the expected value of the left-hand side is 

2 f f e-feW+^K— ) e -^( 2 ) +A )"dnd S 



o jo 



(41) 

(1 - e -fe( 1 )+ A ) f )/( ??j (l) + A) - (1 - e"fa( 2 ) +A ^)/(r ?i (2) + A) 

77, (2) -77, (1) 

Finally, we observe that for every positive integer n, we obtain [recall (29)] 



r]j(n) + X = c z n+ / (1 - (1 -x) n )u z {dx) 
7(o,i] 
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(42) 



c*n + >]( I K- 1 )"" 1 / ^(dx) 

'(0,1] 



and since, r?i 0) (0) = r) z (Q) = 0, rj' z (0) = c z + L G ^xv z (dx) and rj z k \o) 
r x fc ^(dx), for k > 2, it holds that 



fc=0 v ' 



(43) %»+a=>]( ;; )4 fc) (o) 



In particular r]j(l) + A = rj' z (0) = c z + J/ Q ^ x^(dx) and r]j(2) + A = 2^(0) + 
^'(0), so that 77,(2) - % (1) = ^(0) + r£(0). 
To summarize, when EXq = 0, we have 

£X* = 2(i£(0)) a 

x ((1 - er^)M(O) - (1 - e -^W + <(°»')/(2^(0) + rfcO))) 

(44) 

/(v'M+v'M) 

1 _ e -(2^(0)+<(0))t 

"^' (0) 2^(0) + »tf(0) 
which converges to 

2«(0)) 2 - ^(0)<(0) «(0)M(0)) 2 - <(0)/(2^(0)) 



(45) 



^(0)(2^(0) +t#(0)) " 1 + ^(0)7(2^(0)) 



as i — >■ cc. We note that as v z (l,oo) = 0, then clearly whenever either c z > 
or i^(0, 1) 7^ (i.e., Z — N is not identically zero), it holds that 



(46) 



rj' z (0) = c z + [ xu z (dx)> f x 2 v z (dx) = -r ] ' z f (0). 
J(o,i] J [0,1] 



3.5. Levy Z, linear Y and Xq = x. It is of interest to consider the special 
case where Yj = rt for some r > and Xq = x for some x > 0. For the 
case where Z is compound Poisson this model becomes the growth-collapse 
process from [16] where the computation of transient moments turns out to 
be especially tractable. Since 

(47) ^ = X - + t-[ ^dZ s 

r r J(o,t] r 

we may without loss of generality assume that r = 1. Recall (23). Following 
the ideas in the proof of Proposition 3.1 of [4], we first write for a > and 
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b>l, 

b pt / pt \ 6—1 

Ee- aJt ( I e" Js ds) = bEe~ aJt / ( / e~ J " ds ) e~ Ju du 



\Ju 



6-1 

(48) =b I Ee~ a{ - Jt - J ^( j e-^-^cb) e ~( a + b ) J - du 



r-t-U \ 6—1 

b l e -r, 3 (a+b)u Ee -aJ t _ u t I & -J s dfj 



Thus, if T ~ exp(#) for some 9 > and is independent of Z, then since 
the conditional distribution of T — u given T > u is the same as that of T 
(memoryless property), it readily follows that 

(49) Ee~ ajT ( [ T e - Js ds] =— — b — -Ee~ ajT ( F e~ Js ds) . 

For a = we have that, since T\ A T ~ exp(A + 9) and e~ Js 1{ ^ =0 } ds = 

(50) e{[ e^h^ J.)'- (t) | A - e E(£ e-'.l {N ,= 0) d.)"'. 



For a > we have, from the fact that 7\ A T is independent of I^t-^t}) that 

/ /-T \ 6 

i? e - ajT l {WT= o } / e~ J n {Ns=0} t 



o 



(51) 



and thus 



Ee- aJ ^ T l {Tl>T} U^ e- J °ds\ 

a / pT\/\T 



A + 6> 



6 



Ee- aJr l { ^ T=0} / e- Js l {iVs=0} d S 



(52) = £ e "^ 1 {7Vt=0} ^ e~ J ° ds) 



Vj(a + b) + A + (9 
Clearly, when 6 = and o > we have that 

(53) Ee~ ajT l {NT=0} = e ~M«)+ A ) T 



Ee- ajT l {NT=0} { I e~ J °d 



T \ 6-1 







r jj (a) + X + 9' 
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Now 



(54) 



EX% = E 



( xe~ jT 



{N T =0} 



+ 



e~ Ja l 



£ 

k=l 



n 



x k Ee 



-kJ T 



{7V T =0} 



{N a =0} 
T 



ds 



- Js ds 



o 



n—k 



E 



and denoting [recall (43)] 



T 



'1 



{N s =0} 



ds 



(55) fjLi = r )j (i) + \ = c z i+ / (l-(l-x)> 2 (dx) = 2( Mr,f(0), 



'(o,i] ~ \ k 
it follows from (50), (52), (53) and (54), with some manipulations, that 



nl 



EXj> - . 

Ui=i M 



(56) 



1 \k=l 



n 



n 



in 



+n 

i=l 



iH + er 



where an empty product is defined to be 1. Finally, noting that EX^ = 
f™e- 0t dEX? it follows that if {E { \i > 1} are i.i.d. random variables with 
distribution exp(l), then E{/ ^ ~ exp(/Xj). It is well known and easy to check 
that 



( 57 ) n 

i=k 

hence, for 1 < k < n, 



t'i 



e- et dP 



" TP 



(58) n 



i=k+l 



n 

i=k 



t'i 



e~ et dP 



i=k+l 



i=k 



and thus we have the following somewhat curious result. 



Theorem 3. Letpij(t) be the transition matrix function of a pure death 
process D = {Dt\t > 0} with death rates \x%, i > 1 (0 is absorbing). Then 



nl 



(59) 



EX™ - 

n»=i & 



Pn0(i)+E Xfen ^ lM! 



k=l 



Pnk(t) 



nl 



j-in 

Ui=i Mi 



■J5 



1 r 

i * 



Do=n 
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In particular, when x = 0, then 
n ' 

EX? = TfW — p n0 (t) 



(60) = n 



J ■■■ J exp | - ^2 ^i x i j dxi • • ■ dx n 



x-i,...,x„>0 



J" J exp^-t ^2/j,iXij dxi ■ ■ ■ dx n . 



xi,...,a; n >0 



In fact, one may also give a finite simple algorithm with which to compute 
EX™. For the sake of brevity we do it only for the case x = 0. This can be 
done similarly to the Brownian motion in the proof of Theorem 1 on page 31 
of [22] or, equivalently, directly from (60) as follows. Set /o = and for n > 1 
and < a\ < a2 < • ■ • < a n , let 

/„(ai, . . . , a n ) = J ■ ■ ■ J exp f - ^ a;^ J dx x ■ ■ ■ dx n 



xi,...,a;„>0 



i=l 



1 £i=2 x * 

e~ aiX1 dxt 



(61) 



x 2 ,...,a:„>0 







exp - a i x i dx 2 ' ' ' ^ X r 



i=2 



_ /n-l(Q2,- • ■,Qn)-e ai /n-l(«2 ~ «1 , ■ ■ ■ , «n ~ «1 ) 

ai 

Alternatively, if we denote go = 1, and for n > 1 and bi,...,b n > 0, 

(62) ffn(&l, ■ • • > K) = fn{b\M + fe 2 , ■ ■ ■ ,h H h &n)- 

Then 

, co x /, , x 5n-l(fel+fe2,fe3,---,fen)-e" fel S' n _l(fe2,fe3,---,fen 

(63) 5n(0 1 ,...,6„) = - 
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From the above, it is also clear (see also [22], Theorem 1, page 31 for the 
case of a Brownian motion) that, in fact, 

EX^ = t n n\f n (fi 1 t,... } fi n t) 

(64) 

= t n n\g n (fiit, (fi 2 - m)t, . . . , {n n - Hn-i)t) 

is a linear combination of exponentials. An algorithm for computing the coef- 
ficients of this linear combination is equivalent to the above simple algorithm 
which involves only a finite number of additions and multiplications. 

We emphasize that the fact that Theorem 3 holds for all n > 1, and the 
algorithm for the computation of moments, also valid for all n > 1, is special 
for the case where Z is a nonzero subordinator. This is true since this is 
the only case where r]j(n) is finite, strictly positive for all n > 1 and strictly 
increasing. 
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